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Abstract 



We examine existence and stability of relative equilibria of the n- 

vortex problem specialized to the case where N vortices have small 
and equal circulation and one vortex has large circulation. As the 
small circulation tends to zero, the weak vortices tend to a circle cen- 
tered on the strong vortex. A special potential function of this limiting 
problem can be used to characterize orbits and stability. Whenever 
a critical point of this function is nondegenerate, we prove that the 
orbit can be continued via the Implicit Function Theorem, and its 
linear stability is determined by the eigenvalues of the Hessian matrix 
of the potential. For > 3 there are at least three distinct families 
of critical points associated to the limiting problem. Assuming non- 
degeneracy, one of these families continues to a linearly stable class of 
relative equilibria with small and large circulation of the same sign. 
This class becomes unstable as the small circulation passes through 
zero and changes sign. Another family of critical points which is al- 
ways nondegenerate continues to a configuration with small vortices 
arranged in an A?^-gon about the strong central vortex. This class of 
relative equilibria is linearly unstable regardless of the sign of the small 
circulation when A^ > 4. Numerical results suggest that the third fam- 
ily of critical points of the limiting problem also continues to a linearly 
unstable class of solutions of the full problem independent of the sign 
of the small circulation. Thus there is evidence that linearly stable 
relative equilibria exist when the large and small circulation strengths 
are of the same sign, but that no such solutions exist when they have 
opposite signs. The results of this paper are in contrast to those of 
the analogous celestial mechanics problem, for which the A^-gon is the 
only relative equilibrium for N sufficiently large, and is linearly stable 
if and only if A > 7. 
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1 Introduction 



A well known problem in fluid mechanics is that of n point vortices moving 
in the plane, which may be regarded as an approximation to more general 
vortex motion. This is commonly known as the n-vortex problem. Of par- 
ticular interest are solutions that appear fixed when viewed in a uniformly 
rotating or translating frame. These are called relative equilibria and have 
been studied by many authors, see |T2], [1], |9], [13], [15] and [2]. 

The focus of this paper is on relative equilibria of the planar n-vortex 
problem with n = N + 1 where one vortex has large circulation and have 
small, equal circulation. We show that the weak vortices limit to a single cir- 
cle about the strong central vortex as the small circulation tends to zero, and 
that the resulting configuration must be a critical point of a particular poten- 
tial function depending only on the relative positions of the "small" vortices 
on this limiting circle. Conversely, we find that nondegenerate critical points 
of this potential continue to relative equilibria with nonzero circulation, and 
the linear stability of these solutions is determined by the type of critical 
point of the limit potential. 

One case that has been carefully studied has the small vortices in an 
A^-gon around the large vortex in the center. This family is a relative equi- 
librium of the problem regardless of the relative strengths of the central and 
surrounding vortices. Its stability has been studied in detail by Cabral and 
Schmidt in j4j. In particular, they showed that the A^-gon around a central 
vortex is unstable when the circulation of the surrounding vortices is suffi- 
ciently small. A new proof of this is presented here and these techniques 
imply the existence of two more families of critical points of the limit poten- 
tial. One of these other families continues to a linearly stable class of relative 
equilibria of the full problem whenever it is nondegenerate. Members of this 
class of configurations have weak and strong vortices of the same sign. 

The methods used here are parallel to those used to study relative equi- 
libria of the (1 + A^)-body problem with Newtonian inverse square force of 
attraction in celestial mechanics. The results for vortices are, surprisingly, 
very different. For the celestial mechanics problem, the relative equilibrium 
formed by an A^-gon of small masses about a large central mass is stable 
for N > 7 and is the only relative equilibrium for A^ sufficiently large, see 
[I], [H], [in], [E], [S] and [H]. Another difference is that the circulation pa- 
rameter of the vortex problem, which is analogous to mass, can be positive 
or negative. We will see that there can be a dramatic bifurcation as this 
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parameter passes through zero. 

In the next section of the paper we carefully define the problem and build 
the important potential function. In Section |3] we carry out the proof that 
nondegenerate critical points of this potential can be continued to relative 
equilibria with vortices having small, nonzero circulation and one having 
large circulation. Moreover, we show that the linear stability of such a con- 
tinuation is determined by properties of the potential. Section |4] consists of 
applications of Section [3] results to the cases N = 2, 3 and 4. In Section [5] 
we discuss further results for large and do some numerical investigation. 
In particular, we prove that when > 4 the A^-gon consisting of A^ small 
vortices surrounding a strong central vortex is linearly unstable. We use this 
to show that when all critical points are nondegenerate there must exist at 
least two other families of relative equilibria of the full problem. We argue 
that one of these families is made up of linearly stable orbits when the central 
and surrounding vortices are of the same sign. We use numerics to study the 
candidates for the three families and carry out stability calculations. 



2 Relative Equilibria 

Consider an inviscid, incompressible fluid in the plane with vorticity distri- 
bution described by a finite sum of Dirac delta functions. Note that this 
"function" is zero except at the locations of the point vortices. Associated 
to each of these point vortices is an important quantity T{t), called the cir- 
culation of the vortex. Assuming the flow is acted on only by conservative 
forces, we may conclude that T(t) = F, a constant. We therefore treat T as 
a parameter of the problem. 

With this in mind, let qi{t) = {xi(t),yi{t)) G M^, i = I, ■■■,n be the posi- 
tions of n point vortices in the plane with associated circulations Fj. Then 
the equations of motion are described by 

Tiij = ^ (1) 
^ dH 
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where the Hamiltonian function H is 

n 

H{qi,...,qn) = - ^i'^j^^Sh- Qjl- (3) 

For details, see [T^ . 

For n = 2, 3 the system is completely integrable, and all possible relative 
equilibria have been discovered. Proof and examples can be found in [12]. 
Later, we will discuss how these orbits appear in the specialized problem of 
one strong and weak vortices. For n > 3 the problem is nonintegrable, 
and hence more difficult. For proof of the nonintegrability of the four- vortex 
problem, see [16] or [8]. 

We now specialize the general problem to our desired setting. Consider 
iV + 1 vortices go, Qi, ■■■,(1n G with strengths Tq = 1, Ti = e, i = 1, N. 
The equations can be written in the form 

Er.|^. j-o.....iv (4) 

where {x,y)-^ = {—y,x). Next we define what it means to be a relative 
equilibrium of this specific problem and explore some properties of solutions. 

Definition 1. Let {ek}'^=i be a sequence of real numbers such that — >■ 
as A; — 7- oo and let gQ,...,g^ be a sequence of configurations which are 
relative equilibria of the {N + 1) -vortex problem for each k with corresponding 
circulations Fq = l,r^ = Ek, j = 1,...,N. A relative equilibrium of the 
(1 + A^)-vortex problem is a configuration qo,...,q]sf such that q^ — )■ qj as 
A; — )■ oo for each j = 0, A^. 

Remark 1. Both positive and negative circulation e make physical sense. 
The existence and continuation results that follow hold for e of either sign. 
However, linear stability depends on the sign of e, which we point out below. 

Remark 2. To avoid confusion, we remark that the phrase {N + l)-vortex 
problem is meant to refer to the original equations ^ of the n-vortex problem 
with n = A^ + 1. On the other hand, the phrase (1 + A^)-vortex problem refers 
to the limiting problem in which there is one strong vortex and N vortices 
with zero limiting circulation. This is consistent with the celestial mechanics 
literature. 
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Definition [T] does not rule out the possibility that two or more vortices 
collide in the limit. However, we wish to exclude this situation from our 
calculations and so we require that all vortices are bounded away from each 
other by some constant m > 0. We also ignore configurations which become 
unbounded as — > oo, and so we assume that all vortices lie inside a circle of 
radius M about the origin. Both M and m are assumed to be independent 
of k and Sk-, but may depend on A^. Under these assumptions we can prove 
a lemma which asserts that relative equilibria of the (1 + iV)- vortex problem 
have small vortices on a circle of fixed radius about the large vortex, which is 
also fixed. Since the center of vorticity is an invariant of the problem we take 
it to be the origin for simplicity, i.e. we choose go + + ••• + ^n) = 0. This 
will ensure that the large vortex is near the origin. Moreover, by choosing 
an appropriate rotation rate the small vortices will limit to the unit circle. 

The next two results have counterparts in celestial mechanics, see [7j or 

Lemma 1. In a rigidly rotating frame with rotation rate u = 1, all relative 
equilibria, qQ,...,q^, which converge to a relative equilibrium qo,...,qN of the 
(1 + N)-vortex problem satisfy |gg | = 0{ek), Ig^ P — 1 = 0{ek), j = 1, A^. 

Proof. We will drop the superscript k whenever it is unambiguous to do 
so, i.e. when e 7^ is clear. Let q^ = (x^, y^) e M^, = -\- iy'^, and 

^ e^'^^Zj. Suppose ^ = (^q,...,^^) is a fixed point of the system 
written in rotating coordinates. That is, 

^ 1 

= -^u^^o + sJ2-— (5) 
1 ^ 1 

o = -^^e, + 7737^ + ^ 7f ^■ = 1--^- (6) 



For fixed u = 1, 



'^^o = eY,—^ (7) 
i=i ^0 ?i 



N 



1 ^ 
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and so we see that 

N 

iCol < e- (9) 
m 

which proves the first assertion of the lemma. Prom this and the assumption 
that the distances between vortices is bounded below we have 

N 

= 16 - - 61 > 16 - - 161 >m-e- (10) 



m 



Therefore, 



ie.f < |1 - I, + ... + Q + Oie')\ + '^^l (11) 

( MN M(N-1)\ , , 

< 1 + £ ^ + ^ ^ + O £^ 12 

\m — e^ m I 

fMN M(N-1)\ , , 
\ m m J 



Similarly, we find 



fe|2>i_£j^_ ^^(^-l) +0(,2) (14) 

which concludes the proof. □ 

Since any periodic relative equilibrium must rotate rigidly about its center 
of vorticity (assumed to be zero), a necessary condition for a formation to 
be in relative equilibrium is . = 0, j = 1, n. Using this idea, we prove 
the following lemma. 

Lemma 2. Let q = {qo,---,QN) be a relative equilibrium of the (1 + A^)- 
vortex problem and let {r,9) = {ri, ...,rN,Oi, ...9n) be the representation of 
{qi, ...,gAr) in polar coordinates. Then 9 is a critical point of the potential 
function 



^(^) = - (^°^(^^ - + I - 2 cos(^i - 9j))\ 



(15) 



Proof. Let q^ be a sequence of relative equilibria of the full (N + l)-vortex 
problem which converges to as £ — > 0. Since q^ ■ q^ = 0, after suppressing 



the e dependence we have, for j = 1, A^, 

'''' 

( qf + e{qi + ■■■ + gjv)^ A e{qf-qi) \ 
'^''\\q,+e{qr + ...+qN)?^^j^^^ \q, - q^? ) ^ ^ 

eqj ■ {qi + ... + g^v)^ ^ e{qj ■ qi) \ ^^^^ 



\qj + £(gi + ... + g7v)P \qj - qi 



N 

e 



^ '^'"^^(\q,+eiq, + ...+q^)\' \q, - q^^'l ' ^^^^ 



In polar coordinates, an expansion of the above expression in powers of e is 
given by 

^ /I 1 \ 

^eY.r,nsMe, - e.) - ,.^,._,,^,,,„,(,^_,J + (20) 

Due to Lemma [T} dividing both sides by e and taking the hmit as £ — ?■ 
yields 

Therefore, 6' is a critical point of 

V{e) = -Y, (cos(^, - Oj) + ^ log(2 - 2 cos(^, - ^,))'J . (22) 

□ 

Remark 3. In what follows, the Hessian matrix Vqq will play an important 
role. Due to the rotational symmetry of the problem, this matrix will al- 
ways have at least one zero eigenvalue with corresponding eigenvector vq := 
(1,1,...,!)"^. Therefore, as pointed out in [Tl^ . it is appropriate to define 
a critical point, 9, of V to he nondegenerate provided Veg has only one zero 
eigenvalue. Further, we say that 6 is a nondegenerate local minimum (respec- 
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tively, maximum) of the potential V ifVge is positive (negative) semidefinite 
with a one- dimensional null space. 



3 Linear Stability 

Due to the many symmetries and integrals present in our system, the usual 
stability theory for ordinary differential equations does not apply. Here, we 
modify the definitions of nondegeneracy and linear stability to account for 
these degeneracies. 

First of all, it is impossible for a fixed point (or periodic orbit) of the 
system ^ to be hyperbolic because there are always at least four zero eigen- 
values associated to the linearization. Two of these eigenvalues come from 
the invariance of the center of vorticity, which we have assumed is at the 
origin. By defining go '■= ~^(9i + ••• + Qn), the equations for the components 
of go are satisfied once qi,...,qN have been specified. Therefore, we ignore 
these two equations and focus only on the remaining 2N equations, hence re- 
moving two zero eigenvalues from the system. Written in polar coordinates, 
the reduced system becomes, for j = 1, A^, 



d , A f 1 1 



2\ 



— ir^) = e r.rj sin(6'o - OA [ ^ ^ ^ — — | + 0(e , 

dt^ j^. V^l r] + rf -2rjriCos{9j -9i) J ^ ' 

=:eFj{r,e,e) (23) 

^ _ 1 ^ ^ (rj cos(2(g, - 9j)) - n cos(g, - 9^)) _^ ^ 2^ 
dt r| r^{r^ + rf -2rirjCOs{9,-9j)) 

=:Gj{r,9,e). (24) 

The remaining two zero eigenvalues can be readily identified once we linearize 
this system about a relative equilibrium of the problem. Consider a fixed 



point of (23)-(24) which is a member of a sequence of relative equilibria 
converging to a relative equilibrium of the (l+A^)-vortex problem. Then upon 
applying Lemma [T] for e sufficiently small, the matrix M of the linearized 
problem about this fixed point is made up of four N x N blocks: 

/ ^eA + 0{e^) eVee{(t^) + 0{e^) \ 
- I -21 + 0{e) eA + 0{e^) J ^"^^^ 
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where A = (aij) satisfies 

ttij = sin(0j - 0i) j (26) 
flii = ^sin(0j - (27) 

From the form of A and Remark [sj it is easy to see that M{0,vq)^ = 
0{e'^). This computation suggests that one zero eigenvalue may be related 
to the rotational symmetry of the problem. This is true, and we leave the 
details as an exercise. We show in the proof of Theorem [T] that we may ignore 
this eigenvalue by focusing on the complementary eigenspace to the span of 

Vq. 

Due to the Hamiltonian nature of the problem, this eigenvalue must be 
paired with a second zero eigenvalue. Physically, this second eigenvalue is 
born from what can be thought of as the scaling symmetry of the problem: 
associated to any relative equilibrium is a family of relative equilibria which 
are related to the first by a scaling and have different rotation frequencies. 
One may deal with this degeneracy by specifying a rotation frequency, as we 
did in the proof of Lemma [T} 

With this in mind, it is reasonable to define a relative equilibrium of the 
{N + l)-vortex problem to be nondegenerate if the associated linearization 
has exactly four zero eigenvalues (or, once the equations for go are removed 
as described, exactly two zero eigenvalues). 

Aside from the zero eigenvalues, we remind the reader that for Hamilto- 
nian systems asymptotic stability is impossible, since all eigenvalues come in 
pairs with opposite signs. Therefore, one says instead that an equilibrium 
is linearly stable if all eigenvalues are purely imaginary. In the proof of the 
main theorem of this section, we leave the two zero eigenvalues associated 
to scaling and rotational symmetry "as is" and call an equilibrium linearly 
stable if all remaining eigenvalues are purely imaginary and nonzero. 

The first theorem of this section asserts the existence of a convergent 
e-dependent family of relative equilibria whenever there is a nondegenerate 
critical point of the potential. 

Theorem 1. Let cj) = (01, 07v) be a nondegenerate critical point of the 
potential V. Then for p = (1, !,...,!), the configuration {p,(f)) is a relative 
equilibrium of the (1 + N)-vortex problem, that is, there exists a sequence of 
relative equilibria of the {N + 1) -vortex problem which converges to {p,(f)) as 
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e ^0. 



Proof. The proof rehes on two apphcations of the Imphcit Function Theorem. 
Define F(r, ^, e) := {Fi{r, 9, e'h.^^.FUr, 9, e)) and G{r, 9, e) := {Gi{r,9,e), ...,GN{r,9,e)) 
where Fj and Gj are is in (23|)-(p4l. For any G M", set ^ = ^ + 9^,^11 



where 9 G spanlfo}"*" and 6'nuii ^ spanjiio}, Vq as in Remark [sj Then 
F{r,9,e) = F{r,9,e) and DoF{r,9,e) = DoF{r,9,e) because the equations 
depend only on the differences between the angular components. Therefore, 
for the purposes of applying the Implicit Function Theorem, we restrict F 
and DF to spanlfo}"*". 

With this in mind, set = + </)nuii- Since 0) = and the restric- 

tion of DgF{p, (f),0) = Vee(0) to this subspace is invertible, we obtain open 
sets U, V and W containing r = p, 9 = (p and e = 0, respectively, and a con- 
tinuously differentiable function / : W x W — )■ V such that F(r, /(r, e), e) = 
whenever (r, e) G W x W. This implies that there exists a sequence of so- 
lutions of the (A^ + l)-vortex problem which converges to (p, 0) as e — > 0. 
However, at this point it is not clear that these solutions are relative equi- 
libria of the problem, since it is possible that different vortices have different 
rotation rates. 

To finish the proof, we apply the Implicit Function Theorem to G. Set 
9 = f{r,e) with (r, e) restricted to U x W, as determined in the previous 
paragraph. Then G{p, f{p, 0), 0) = G{p, 0, 0) = 1 and 

A.(G(r,/(r,£),£))|(^_^^^(^_^^ = D,G(p,0,O) + |^(p, 0)D,G(p, 0, 0) (28) 

= -2L (29) 

Thus there exist open sets Ui C U, Wi C W containing r = p and e = 0, 
and a continuous function g : Wi — > Ui such that G{g{e), f{g{e),e),e) = 1 
whenever e G Wi. It is also true that F{g{e), f{g{e), e), e) = for e G Wi 
and so the proof is complete. □ 

The symplectic structure associated to any Hamiltonian system provides 
a rich theory and many tools for analysis. Before stating the main theorem 
of this section, we list one definition and result that will be used in the proof. 
However, we will not go into the details of the theory here. For further 
discussion and proof of the lemma, see pT] . 
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Definition 2. Let J be the 2N x 2N block matrix given by 



The skew inner product of two vectors defined to be 

n{v, w) = v^Jw. (30) 

Lemma 3. Suppose each eigenvector v of M satisfies the inequality Q{v,v*) ^ 
0. Then all eigenvalues of M are purely imaginary. 

Theorem 2. Suppose {p'^,(f)'^) = (pf , p|^, 0f , 0^) is a sequence of rela- 
tive equilibria of the {N + l)-vortex problem which converges to a relative 
equilibrium {p,4>) = (1, 1, 0i, 0Ar) of the (1 + N)-vortex problem as 
e — )■ 0, where (p is a nondegenerate critical point of the potential V defined by 



Equation (20). Then {p'^,(f)^) is nondegenerate for e sufficiently small. 
Moreover, in this case the configuration (p^, 0^) for e > is linearly stable 
if and only if (p is a local minimum of V. Likewise, for e < it is linearly 
stable if and only if (j) is a local maximum of V . 

Proof. Much of the foUowing proof is modeled on analogous arguments made 
by Moeckel in [TU] for the corresponding celestial mechanics problem. Let M 



denote the linearization matrix (24) about the relative equilibrium [p'^^cjf). 
Throughout the proof we suppress the explicit e-dependence when there is 
no ambiguity. Let A G C be nonzero and consider the matrix M — \I2Nx2N'- 



-21 + 0{e) -\I + eA + 0{e^) 

Since the matrix Vgg{(f)) is assumed to be nondegenerate, it has exactly one 
zero eigenvalue. This, coupled with the discussion prior to Theorem [T] implies 
that M has at least two zero eigenvalues. Choose < ci < 1 small enough 
that Id > ci for any nonzero eigenvalue, (, of Vgg^cj)). We now consider 
det(M — A/) restricted to the region ciy/e < |A| < e^^^. An application of 
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Lemma |4] yields 



det(M - A/) = det(-A/ + 0{e^)) det(-A/ + 0{e) 

+ {-21 + 0{e)){-XI - 0{e))-\6Vge + 0{e'))) 

= det(-A/ + 0{e)) det(-A/ - '^Vee + 0{e)) (31) 

A 

= det(/ + C(v^)) det(A2/ + 2eVee + 0{e^^*)) 
= (1 + 0(v^)) det(A'/ + 2eVee + 0(6^^^)). 

The only way that the right side can be zero in this region is if A = 0{^/6). 
Set A(e) = y/e'y{e) and define lim£_j.o7(£) =: 7o- Then 

det(£72j + 2eVee + 0(6^/^)) = det(72/ + 2Vee + 0(v^)). (32) 

Now observe that det(M — XI) is a polynomial of degree 2N in A, call it 
Pe(A). Moreover, 

lim ^PeiV^i) = det(7o'/ + 2Vee). (33) 

The latter is a polynomial of degree in ^ = 7q which has — 1 nonzero 
roots. Therefore, for e sufficiently small there are 2N — 2 zeroes of of the 
form A = y/e'y{e) such that •y^e) ±70 and \j{e)\ > Ci. As we remarked 
above, M has two zero eigenvalues, so those two, together with the 2N — 2 
just found comprise the entire spectrum of M. 

The preceding arguments hold for e ^ 0, regardless of sign. For the rest of 
the proof we assume e > and remark that completely analogous arguments 
can be made for e < 0. 

Suppose that e > and is not a local minimum of the potential. Then 
2Vg0{(j)) must have a negative eigenvalue which implies that for all e suf- 
ficiently small, 7(e) must have nonzero real part. It follows that M has 
eigenvalues with nonzero real part for e small and so the relative equilibrium 
is not linearly stable. 

The final piece of the proof is to show that in the case that e > and </> 
is a local minimum of the potential, the hypothesis of Lemma [3] holds. Let 
V = {vr, ve)'^ be an eigenvector of M associated to the eigenvalue A = y/e'y{e). 
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Using the equation (M — \I)v = we see that 

- Xvr = 0{e) (34) 

which imphes that Vr = 0{y/e). Set Vr = y/ewr{e) so that the second 
component of the eigenvector equation becomes 

-2y/iwr-Veive = 0{e). (35) 

From this we see that Wr = —^vg + 0{^/6). Observe also that if we impose the 
normahzation condition |frP + |weP = 1 we can conclude that l^gp = l+0{e). 

Suppose that C G is a nonzero eigenvalue of 2Vg0. It follows from the 
assumption that is a local minimum that ^ > 0, and so —7^ = C + o(l) 
then implies that 7 = /3z + o(l) where /3 G M is nonzero. 

Finally, a short computation gives 

fl{v, v) = 2lm{vr ■ ve) 

= 2Im(^^hp + 0(£)^ (36) 
= P^e + 0{e) 

where the last equality follows from the normalization assumption. There- 
fore, for e small and positive the hypothesis of Lemma [3] is satisfied and the 
desired result follows. □ 

Remark 4. According to the theorem, relative equilibria which are near non- 
degenerate local maxima or minima of the potential undergo a complete loss 
or gain of stability as the circulation strength e passes through zero, while so- 
lutions near saddle points remain unstable. Thus, the behavior of the system 
if very different depending on whether or not the circulation of the "small" 
vortices has the same sign, or opposite sign, as that of the large central vortex. 

4 Results for Small N 

In this section we describe solutions of the (l+A^)-vortex problem for = 2, 3 
and 4. For N = 2 and 3, we discover all possible relative equilibria of the 
problem and examine linear stability. For iV = 4 we use a numerical root 
finder based on Newton's method to find critical points of the potential and 
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also examine stability. 
4.1 N = 2 

It is well known that the only relative equilibria of the 3-vortex problem are 
collinear and equilateral triangle configurations. First consider the triangle. 
We compute the Hessian matrix of V{9i, 62) with, e.g., 61 = 0, 6*2 = |. Then 

^.(0, i) - ( j| 7 ) 

which has eigenvalues 3 and 0. Therefore, for e > Theorems [T] and |2] imply 
the existence of a linearly stable family of relative equilibria converging to this 
configuration as e — 0. Similarly, we discover that the collinear configuration 
corresponds to a linearly unstable family. 





(a) (b) 

Figure 1: Relative equilibria of the 1+2 vortex problem. For e > 0, configu- 
ration (a) is the limit of a family of linearly unstable relative equilibria and 
(b) is the limiting configuration of a sequence of linearly stable equilibria. 
For e < 0, the two stability types switch. 



4.2 iV = 3 

By straightforward computations using the potential one can show that the 
(1 + 3)- vortex problem has exactly three relative equilibria up to rigid rota- 
tions and permutations of indices. Figure [2] illustrates these three families. 
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37r 
4 



In Figure 2a, the vortices on the left side of the circle each form an angle of 
with the positive x-axis. This configuration is a saddle point of the potential 
and so it corresponds to a linearly unstable family of solutions for both e > 
and 6 < 0. The equilateral configuration about the central vortex is a local 



maximum of the potential. In Figure 2c , the small vortices are separated by 
an angle of 7r/4 and this configuration is a local minimum of the potential. 



(a) 



(b) 



(c) 



Figure 2: Relative equilibria of the 1+3 vortex problem. For e > the 
configurations (a) and (b) are limits of a family of linearly unstable orbits, 
while configuration (c) corresponds to a family of linearly stable orbits. For 
e < 0, configuration (b) is the limit of a family of linearly stable relative 
equilibria, while the other two correspond to unstable solutions. 



4.3 iV = 4 

The graphics in Figure |3] show the numerically observed relative equilibria 



for = 4. In configuration 3a, two of the vortices form an angle of | with 
the negative x-axis. The remaining two small vortices are at 1 and vr. Using 
the potential, we calculate the nonzero eigenvalues of the Hessian at this 
configuration to be 4, —|, 1 and so this is a saddle point. We next calculate 
the nonzero eigenvalues associated to the A^-gon to be 2, — |, — |. Thus the 
A^-gon is also a saddle point of the potential. In fact, we will show in the next 
section that the A^-gon is always a saddle point of the potential for A^ > 4. 
The vortices in the third configuration are separated by an angle slightly less 
than |. The eigenvalues are approximately 12.4,8.4, and 3.7 and so we find 
this to be a local minimum of the potential. 
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(a) (b) (c) 

Figure 3: Numerically observed relative equilibria of the (1+4) -vortex prob- 
lem. For e > 0, they correspond to an (a) unstable, (b) unstable, (c) stable 
family of relative equilibria. For e < 0, all three of the converging families 
are linearly unstable. 

We would like to emphasize the fact that all three families of relative 
equilibria become unstable as e passes through zero and becomes negative. 
This result also applies for larger values of N, which we describe in the next 
section. 

5 Results for Large N 

In this section we prove that the N-gon is always a nondegenerate saddle 
point of the limit potential if > 4. This implies that the limit potential has 
additional critical points, one of which must be a minimum. In the simplest 
case, when all critical points are nondegenerate, there must be at least three 
distinct families. In numerical searches, we find just these three families 
for > 4, and they are nondegenerate in every case computed. When 
nondegenerate. Theorem 2 implies these critical points continue to e 7^ 0. 
The minima continue to linearly stable equilibria for e > 0. When e < 0, 
the relative equilibria corresponding to both the A^-gon and the minimum 
become highly unstable. We investigate the hypothesis that this is also the 
case for the third observed family, thus highlighting the dramatic loss of 
stability that occurs when the circulation of the "small" vortices is opposite 
in sign from that of the "large" central vortex. 
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5.1 Instability of the N-gon and Implications 

For all N, the A'"-gon about a central vortex at the origin is a relative equi- 
librium of the (1 + A'^)-vortex problem. In the previous section we explained 
that when iV = 2 or 3 any sequence of relative equilibria converging to the 
iV-gon is made up of linearly unstable orbits for 5 > and linearly stable 
orbits if £ < 0. When = 4 the iV-gon is a saddle point of the potential and 
so for both £ > and £ < it corresponds to families of hnearly unstable 
periodic orbits. We now show that this last statement generahzes to > 4. 

Theorem 3. For N > 4, the Hessian matrix of the potential V{9) eval- 
uated at the N—gon has at least one positive eigenvalue and one negative 
eigenvalue. 

Proof. The Hessian has the form of a circulant matrix with entries 

V,,, = -cos(«.-«,)-^-^J^^. (37) 
V,,,. = -J2^e.e,- (38) 

Let (ao,ai, ...,aAr_i) denote the first row of this matrix. The eigenvalues of 
such a matrix can be computed by evaluating the generator polynomial 

q{t) = ao + ait + C2t^ + ... + a^-it^'^ (39) 

of the matrix t — , j = ...,N — 1 and u — exp ^ is the primitive A'^th 
root of unity. Observe first that g(l) = 0. Next we evaluate q{u): 

q{iu) = -J2 Ve,e, + ujVe,e, + co^Ve.es + - + co'^-'Ve.e^ (40) 

= Vg,e,ico - 1) + Ve,e,iu;^ - 1) + ... + Vg,0^{uj^-' - 1) (41) 
^J2^e,o,{uj^-'-l) (42) 

^y(- cos(^, - e,) ^— -"l fexp ^Al^^ _ . 

^ V ^ 2-2cos(^i-^,); V ^ J 

(43) 

In the case of the N-gon, if we assume that the angles are ordered counter- 
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clockwise around the unit circle with 9i = then 9j — 9i = ^^"^jp^- Further, 
since this matrix is Hermitian all of its eigenvalues must be real and so we 
can ignore the imaginary part of the sum. Therefore, we obtain 

Finally, we compute ^(co'^) which gives 

whenever > 4. □ 

The theorem was proved by computing the first three eigenvalues. How- 
ever, it will be useful for what follows to produce estimates on more eigen- 
values. We do this now. In general, for < j < — 1, 

(47) 

We examine each of the three sums separately. Observe first that 



fe=l 



As in the proof of Theorem |3] we ignore the imaginary parts of the sums. 
Next observe that the third term satisfies 
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for j = 2, N — 2. Now write the second term as follows: 

k=l ^ ^ fc=l 

We can use geometric series to compute this sum, except when j = 1 , — 1 
and when j = y + 1, y — 1 if is even. For j = 1 and — 1, the entire 
sum is equal to — |. For all other j, we compute the second term, above, to 
be —1, even in the special cases listed. Therefore, we find that Vge at the 
A^-gon has one zero eigenvalue, two eigenvalues which are — | and all other 
eigenvalues, A = q{uj^), satisfy 

q{con = -1 - (-1) + bU, N) = b{j, N) > 0. (51) 

Note that since the A^-gon has only two negative eigenvalues, its continuation 
has only two directions of instability when e > 0. However, when e < these 
correspond to two directions of stability, while the A^ — 3 positive eigenvalues 
correspond to directions of instability. Thus, for large A^ solutions near the 
A^-gon are much more unstable when 5 < 0. 

We now compare Theorem [3] to the results in In the article, the 
authors showed that the A^-gon about a central vortex is linearly, and in fact 
Liapunov, stable for a bounded interval of possible strengths of the central 
vortex and that it is unstable whenever the strength is outside of this interval. 
In particular, they showed that if F is the common strength of the vortices 
lying on the A^-gon, and pT is the strength of the central vortex, then the 
configuration is Liapunov stable if and only if A^ satisfies 



N^-8N + 8 (AT -1)2 

— <p < - — - — — for N even, (52) 

+ ^ <p<^j^forNodd. (63) 

16 4 

In the setting of this paper, F = e and so p = -. Thus for e sufficiently 
small (as guaranteed by the Implicit Function Theorem in our arguments), p 
does not fall into the interval for stability and so the configuration is unstable. 
This agrees with our result. Moreover, it is interesting to note that their 
result forces an upper bound of O {-^) on e. See ?? for related results in 
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the (1 + A^)-body problem. 

We finish this section with an important corollary to Theorem |3] 

Corollary 1. Assuming all critical points of V are nondegenerate, there 
are at least three distinct families of relative equilibria of the (1 + N)-vortex 
problem for N > 3, one of which continues to linearly stable equilibria when 
e>0. 



Proof. The case = 3 was considered in Section For > 4, we showed 
that the A^-gon about a central vortex makes up one of these families, and 
that it is linearly unstable for e sufficiently small. We now show that 
V has a local minimum. This family continues to linearly stable relative 
equilibria when £ > by Theorem |2} Without loss of generality, assume 
^1 = and < ^2 < ^3 < ••• < ^jv < 27r. All other critical points of the 
potential are rotations or permutations of critical points of this form. Then 
is a continuous function on the interior of the specified "wedge" of the 
A^-cube (0, 27r)^ that tends to positive infinity at the boundaries (e.g. as 
02 — >■ ^3)- Thus V must achieve a local minimum inside the wedge. 
With restrictions on 6 as above, consider the transformation 

(^^1, 02, ...,9n) ^ (^1,^2 - Ou ...,9n- 9n-i) =■■ ivu V2, Vn) (54) 

where r/i = 0, rji > 0, i = 2, N and ?7i + ^72 + ••• + '7Ar < 27r. Since V{ri) does 
not depend explicitly on rji we simply drop this coordinate and consider the 
function of A^ — 2 variables V{7]2, rj^-i)- This removes the zero eigenvalue 
from the associated Hessian. Under this transformation, the signs of the 
remaining eigenvalues of the Hessian are preserved. Then in the case of the 
A^-gon, Vrjr^ has exactly 2 negative eigenvalues. Therefore, the index of the 
A^-gon for the gradient vector field is one. The index of the minimum is 
also one. Thus we may apply the Hopf Index Theorem (see, for instance, 
the standard text j3]) to and obtain a third critical point of V{ri) having 
index negative one. This implies the existence of a third critical point in 
the original coordinates which is a relative equilibrium of the (1 + A^) -vortex 
problem and is distinct from the A^-gon and the minimum. □ 

We remark that if the minimum is isolated then it can be continued by 
topological arguments to the e 7^ case. While more degenerate situations 
are conceivable, all the critical points seen numerically have been nondegen- 
erate. Assuming this is true of the third family given by Corollary [T| note 
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that the associated Hessian matrix must have an odd number of negative 
eigenvalues, and so it must be a saddle point of the potential when is odd. 
In this case, Theorem |2] implies that its continuation to a relative equilib- 
rium of the full problem is linearly unstable for all e 7^ sufficiently small. 
Numerics suggest that for all N, even or odd, this third critical point is non- 
degenerate with exactly one negative eigenvalue. For this reason we expect 
its continuation to be unstable for e > and highly unstable for e < 0. 
Thus, there is evidence that all three families of relative equilibria become 
very unstable as e becomes negative. 



5.2 Numerical Observations 

We supplement the discussion in the previous subsection with numerical find- 
ings. Using an algorithm based on Newton's method, we were able to locate 
three distinct families of critical points of the potential. We have run this 
algorithm for < 100, and in all of these cases the critical points are relative 
equilibria of the (1 + A^)-vortex problem. In other words, they continue to 
relative equilibria of the full problem with e ^ sufficiently small. 

A representative from each of the three families is shown in Figure |4] for 
A^ = 25. Compare these equilibrium types to those in Figures [1} |2] and |3] 
It is natural to wonder if these three critical points correspond to relative 
equilibria of the full problem for all A^, and if they are exactly those critical 
points guaranteed by Corollary [T] We cannot answer this question conclu- 
sively at this time, but we have numerically observed only these three families 
for A^ < 100. 

For the two relative equilibria with "clusters" of small vortices, we have 
observed that as A^ increases the cluster tends to fill out the unit circle more 
and more, and appears to approach the A^-gon. This filling in process is 
beginning to become apparent in Figure |4j 
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(a) (b) (c) 

Figure 4: Three relative equilibria of the (1 + 25)-vortex problem. 

Using a modification of the algorithm, we were able to identify one of the 
three classes as containing nondegenerate minima of the potential. Therefore 
this class corresponds to linearly stable sequences of relative equilibria of the 
full problem when e > 0. Three representatives are shown in Figure |5} 




(a) (b) (c) 

Figure 5: Some relative equilibria of the (1 + A^)-vortex problem which are 
minima of the potential, (a) N = 30, (b) = 45, (c) = 60. 



To investigate the stability of the third class of critical points we computed 
eigenvalues of the associated Hessian matrix for N < 100. We discovered that 
the matrix has exactly one negative eigenvalue and N—2 positive eigenvalues, 
thus leading us to believe that, like the A^-gon, nearby relative equilibria are 
highly unstable when e < 0. We reiterate that this third family of relative 
equilibria has an odd number of negative eigenvalues (at least for < 100) 
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and is therefore a candidate for the third family of critical points guaranteed 
by the Hopf Index Theorem in the proof of Corollary [Tj 

Finally, we remark that for all three observed families of critical points, 
the positive eigenvalues of the Hessian appear to increase monotonically with 
N. This is still more evidence for the striking instability of all three relative 
equilibria of the full problem when e < 0. 

6 Conclusions 

In this article, we examined solutions of a specialized problem in point vortex 
dynamics: the (1 + iV)-vortex problem. We were able to show that all solu- 
tions of this problem must be critical points of a specified potential function 
and that whenever such a critical point is nondegenerate in some appropri- 
ate sense, there exists a nondegenerate family of relative equilibria of the full 
(A^ -|- l)-vortex problem which converges to it. Moreover, we were able to 
exploit properties of Hamiltonian systems to show that the linear stability of 
a member of such a family is intimately tied to the potential. 

We used these results to prove that the A^-gon about a central vortex 
is a linearly unstable relative equilibrium when the "weak" vortex circula- 
tion is sufficiently small if > 4. We further showed that when all critical 
points of the limit potential are nondegenerate, there are at least two other 
distinct families of relative equilibria with vortices having small, nonzero 
circulation. One of these families was shown to be linearly stable when the 
central and surrounding vortices are of the same sign. Our further numerical 
investigation supports the hypothesis that the problem has no linearly stable 
relative equilibria when the weak and strong circulations are opposite in sign. 

Acknowledgement: The research of A.M.B. was supported by the Center 
for BioDynamics at Boston University and the NSF (DMS0602204 EMSW21- 
RTG). C.E.W. was supported in part by DMS-0908093. He would also like 
to thank Paul Newton for a very useful discussion of relative equilibria of 
point vortices. 

A Appendix 

The following lemma was used in the proof of Theorem [2j 
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Lemma 4. Let A, B, C , and D he real n x n matrices and define 

M:=(^ d)- (55) 

Then detM = dct(A) det{D-CA-^B) = det{D) det{A- BD'^C) whenever 
A and D are invertible. 

Proof. The result follows immediately from the following observation: 

n_^;f-.J (56) 

(57) 

□ 



C I J \ D- CA-^B 

I B \ f A- BD-^C 
D [ D-^C I 
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